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i— |" Abstract 

43 : 

^ . In this work we present a simple thick braneworld model that is generated by 

Q-ij an intriguing interplay between a 5D cosmological constant with a de Sitter metric 

induced in the 3-brane without the inclusion of scalar fields. We show that 4D 
gravity is localized on this brane, provide analytic expressions for the massive Kaluza- 
! Klein (KK) fluctuation modes and also show that the spectrum of metric excitations 

J> ' displays a mass gap. We finally present the corrections to Newton's law due to these 

massive modes. This model has no naked singularities along the fifth dimension 

00 . 

\q , despite the existence of a mass gap in the graviton spectrum as it happens in thick 

branes with 4D Poincare symmetry, providing a simple model with very good features: 

Q\ ■ the curvature is completely smooth along the fifth dimension, it localizes 4D gravity 

and the spectrum of gravity fluctuations presents a mass gap, a fact that rules out 
| the existence of phenomenologically dangerous ultralight KK excitations in the model. 

We finally present our solution as a limit of scalar thick branes. 

•rH . 

X: 

The model and its physical consistency 

Within the framework of the braneworld models embedded in a spacetime with extra di- 
mensions and after the success of the thin brane models, where singularities are present at 
the position of the branes [l]-[6], in solving the mass hierarchy and 4D gravity localization 
problems, it has became a matter of interest to find smooth braneworld solutions. In some 
models, such solutions are obtained by introducing one or several scalar fields in the bulk 
[7]-[17] and the large variety of scalar fields that can be used to generate these models gives 
rise to different scenarios without singularities at the position of the branes. By following 
the line of smoothing out the brane configurations we show that 4D gravity can be localized 
on a de Sitter 3-brane generated in five dimensions by gravity with a positive cosmological 
constant only, avoiding at all the use of scalar fields. Thus, our thick brane is not made 
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of scalar matter but rather is modeled by an intriguing relation between the curvatures 
generated by the 5D and 4D cosmological constants. 

In this scenario, the inclusion of the 5D cosmological constant renders a completely 
smooth 5D manifold in the presence of a mass gap in the KK graviton spectrum of metric 
fluctuations, in contrast with previous results found in the context of thick braneworlds 
with 4D Poincare symmetry [8], [14]. 

The complete action for the braneworld model is expressed as follows 



S = 2M 3 J d 5 a:^ (R - 2A 5 ) . 



(1) 



It describes 5D gravity with a bulk cosmological constant A 5 , and M is the gravitational 
coupling constant in 5D (in fact we can write ^ instead of 2M 3 ). 

The Einstein equations with a cosmological constant in five dimensions are given by 



Gab = —A-59AB 



(2) 



For the background metric we use the ansatz of a warped 5D line element with an induced 
3-brane with spatially flat cosmological background that reads 



ds 2 



,2/(<r) 



-dt 2 + a 2 (t) (dx 2 + dy 2 + dz 2 ^j + da 



(3) 



where f(a) is the warp factor and a(t) is the scale factor of the brane. 

By using the ansatz (3) we can compute the components of the Einstein tensor 
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where " /" and " • " are derivatives with respect to the extra dimension and time, respectively. 
The set of equations (2) can be rewritten in a simple way: 
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f' 2 = 
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(5) 
(6) 



If we derive the restriction equation (6) with respect to the extra dimension in order 
to compare it to (5), we get a second order differential equation for the scale factor with a 
general solution of the following type 



a(t) 



ce 



Ht 
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where c and H are arbitrary constants. Thus, mathematical consistency of the Einstein 
equations forces us to have a scale factor corresponding to a de Sitter 4D cosmological 
background defined by a(t) = e Ht , since the constant c can be absorbed into a coordinate 
redefinition. Once we have determined the form of the scale factor, the solution for the 
warp factor is straightforward and reads 

/((7) = ln(yC0s(k7)), (7) 

where — n/2 < a < n/2, the constant b is inversely proportional to the thickness of the 
3-brane and is determined by the 5D cosmological constant as follows: 

6 2 = (8) 
By performing the following coordinate transformation 

z = J e- f(a) da, (9) 
the 5D metric (3) gets an overall conformal warp factor 

ds 2 = e 2f(z) [g^dx^dx" + dz 2 ] , (10) 

or, in terms of the metric qmn = ^ 2 ^9mn, which enables us to easily compute the 5D 
curvature scalar R in terms of the R and the warp factor through the known formula [18] 

R = e~ 2S [R - 2(d - l)g MN VMVNf - (d - 2)(d - l^WW] , (11) 

where d is the total number of dimensions. This fact, in turn, allows us to separate an 
effective 4D action from the 5D one 

S eff D J d 4 x^g~ |2M 3 J°° [e 3f R + 4H 2 e 3f (5 sech 2 (^) - 3) - 2e 5/ A 5 ] dz} . (12) 

After integrating over the extra coordinate z(a) we shall make a direct comparison of the 
result to the 4D Einstein-Hilbert action on the brane 

S brane = 2M 2 pl J d A x^g ( 4 i? - 2A 4 ) . (13) 

Thus, in order to derive the scale of 4D gravitational interactions we look at the curva- 
ture term and perform the following integration: 



M 2 Pl 



= M e^dz = — seAHHzyj, = = ^ , (14) 

which turns out to be finite, as it should be for a well defined 4D theory. 
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The second and third terms in the integral (12) contribute to the definition of the 4D 
cosmological constant A 4 



4M 2 A 4 



8M 3 H 2 e 3f (3 - 5sech 2 (Hz)) + 4M 3 e 5/ A 5 ] dz. (15) 

3 L \ / J 

After integrating along the fifth dimension and simplifying the resulting expression we get 

A 4 = 3H 2 . (16) 

This relation can also be obtained from the Friedmann equations derived by directly varying 
the action (13). 



Metric perturbations, localization of 4D gravity and 
corrections to Newton's law 

In order to compute the corrections to Newton's law we need to perform some preliminary 
work. First of all we perturb the ansatz for the metric and then we compute the Einstein's 
equations to first order. For the sake of simplifying the calculations we use the freedom we 
have to impose an axial gauge h^M = on the perturbed metric: 

ds 2 = e 2fia) [g^ + V] dx"dx u + da 2 . (17) 

After imposing the transverse traceless condition = d^h^ = 0, the system of per- 
turbed Einstein equations to first order adopts the form: 

dl + 4/0, + e- 2 ^ (-d 2 - 3Hd t + e- 2Ht V 2 - 2H 2 )] h, v = 0, (18) 

where V 2 is the Laplacian in flat space and h^ u are the transverse traceless modes of the 
metric fluctuations. 

In order to simplify calculations we perform the coordinate transformation (9) which 
allows us to study in detail the physics of the previous equation. Thus, the equation (18) 
takes the following form 

'd 2 + 3f'd z - d 2 - 3Hd t + e~ 2Ht V 2 - 2H 2 ] = 0. (19) 

At this point we shall make use of the following separation of variables 

V = e'^ z ^^(z)g(x). (20) 



function ^^(z) have been omitted for convenience) for the 5D sector 



This change of variables leads to a Schrodinger-like equation (where the indices of the 

for conve 

(-d 2 + 9 -f' 2 + 3 -f"-m 2 )nz) = 0, (21) 
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where one can identify the quantum mechanical potential Vqm as follows 

V QM = \f 2 + \f. (22) 

Finally, the 4D equation that we obtain reads 

{-d 2 - 3Hd t + e~ 2Ht V 2 - 2H 2 ) g(x) = -m 2 g(x), (23) 

where the integration constant m 2 is defined as the mass in a 4D de Sitter spacetime [19]. 
The warp factor in the language of the variable z is given by 

/ = In ^jsech(Hz)^ . (24) 

Thus, in terms of this coordinate, our quantum mechanical potential takes the form of 
a modified Poschl-Teller one 

V QM = -^H 2 sech 2 (Hz) + ^H 2 , (25) 

which ensures the existence of a mass gap determined by its asymptotic value, namely, by 
jH 2 , or equivalently m = ^j-. 

Thus, we have a classical eigenvalue problem for the following Schrodinger equation 

(-d 2 z - ^-H 2 sech 2 (Hz) + ^H 2 - m 2 ^j V(z) = 0, (26) 

with the general solution given by 

(*) = Ci P% (tanh(#z)) + C 2 Q% (tanh(^)) (27) 

2 2 

where P% and Q^i are associated Legendre functions of first and second kind, respectively, 



of degree v — 3/2 and order \x = y | — ^ [12]. 

Thus, the last differential equation has two discrete bound states: the first one corre- 
sponding to a massless bound state with energy E = —\H 2 and /i = 3/2 

y (z) = k sech%(Hz), (28) 

physically interpreted as a stable graviton localized on the brane, and a second one corre- 
sponding to an excited state with m 2 = 2H 2 , energy E\ = —\H 2 and /i — | 

tfi(z) = k 1 wnh(Hz)aeeh*(Hz), (29) 

which represents a massive graviton also localized on the brane; finally we have a tower of 
massive continuous modes that start from m > ^j- described by the following eigenfunctions 
with imaginary order ji — ip [12] 



tf m (z) = C x Pf lp (tanh(tfz)) + C 2 Qf p (tanh(tfz)) , (30) 

2 2 
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where p = ~ §, that must asymptote to plane waves. 

This result prepares to us the way to compute corrections to Newton's law, since the 
massive modes contribute to Newton's law with small corrections coming from the fifth 
dimension. In order to find these corrections we consider 2m > 3H and C 2 = in equation 
(30) since the continuous spectrum of eigenfunctions asymptotically describes plane waves 

= C±Pf ip (tanh(ifz)) ~ -^e ±tHpz . (31) 

2 V27T 



After this, we further take the thin brane limit H i — > oo and locate a probe mass in the 
center of the brane in the transverse direction. The corrections to the potential generated 
by massive gravitons can be expressed as follows [9] 

U( r ) ~ (G 4 + M- 3 e" mir |^i(^ )| 2 + M~ 3 f°° dme~ mr V (m \z ) (32) 

T V J mo / 

= MlMl ( G4 + AG 4 ) , (33) 
r 

where the brane is located at z = zq, G^ is the gravitational coupling constant in 4D, 
^i{zq) is the wave function of the first excited state and since it is an odd function, it does 
not contribute in the thin brane limit that we are considering, finally, ^/^(zq) represent 
the continuous massive states that need to be integrated over their masses. Thus, the 
corrections to Newton's law contributed by the massive continuous states have the following 
form (see [12] for details): 



AG a ~ m; 



e-* Hr (. „ / 1 



r(-M) 2 



1 "(77;)) < 3J ) 



indicating its small character. 

It is worth mentioning that in our setup there also exists a relation between the ground 
bound state ^§{z) = k,Qez^ z \ which ensures 4D gravity localization, the smooth character 
of the 5D Ricci scalar R and the behaviour of the quantum mechanical potential Vqm in 
the same spirit as it was done for flat thick 3-branes in [14], where it was shown that if 
one requires the existence of a mass gap in the KK graviton spectrum along with the local- 
ization of 4D gravity, then the 5D manifold necessarily develops naked singularities at its 
boundaries. Conversely, if one imposes a smooth character of the Ricci scalar together with 
4D gravity localization, then the quantum mechanical potential asymptotically vanishes, a 
fact that rules out the presence of a mass gap in the graviton spectrum of KK excitations. 
In this work we include a positive bulk cosmological constant A 5 together with an induced 
flat FRW metric (of de Sitter type, in particular) on the 3-brane. 

It is easy to show that for the setup (1) with the ansatz (17) and for the potential Vqm 
considered in (25), for smooth configurations at the position of the branes we can write the 
following relation 




(35) 
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which apart from the quantities included in [14], it involves the scale factor. Thus, from this 
expression we see that if one requires the presence of a mass gap in the graviton spectrum 
of KK excitations, the potential Vqm must asymptotically adopt a positive value; if indeed 
one requires 4D gravity to be localized, which means that — > as z — > oo, then the 
Ricci scalar can be either finite or singular depending on the asymptotic behaviour of the 
expression in brackets: if it asymptotically vanishes, then the Ricci scalar can be finite along 
the fifth dimension, while if it asymptotically tends to a finite or infinite value, then the 
curvature scalar will possess naked singularities at the boundaries of the manifold. A similar 
result was recently obtained in [20] for AdS5 bulk geometries with both dS and AdS induced 
metrics on the 3-brane, where some examples are displayed as clarifying illustrations. 

In our case the scale factor has the particular form corresponding to a de Sitter metric 
a = e Ht , yielding 

i 

(36) 



2 4 

3H - -Vqm 



since the quantum mechanical potential Vqm asymptotically approaches the positive value 
9H 2 /4, ensuring the presence of a mass gap in the graviton spectrum of KK fluctuations, 
then the expression in brackets asymptotically vanishes, implying that the scalar curvature 
is regular along the fifth dimension, which is indeed the case since it is constant R = 10A 5 /3. 

Thus, this result guarantees the smoothness of our braneworld obtained without the use 
of scalar fields even in the presence of a mass gap in the spectrum of KK graviton fluctu- 
ations, an important phenomenological fact that rules out the existence of the dangerous 
arbitrarily light KK excitations. 



The solution as a limit of scalar thick branes 

The above presented solution (7) can also be obtained as a special case or limit of scalar field 
thick branes, with a bent de Sitter 4d metric, by fixing one parameter of these solutions. 

By considering the action of a bulk scalar field with a self-interacting potential V(<p) 
minimally coupled to 5d gravity 



d 5 xJ\G 



\r 5 



i(V0) 2 -V(0) 



(37) 



and assuming a de Sitter 3-brane in the metric ansatz (3), one gets the following solution 
reported by Gremm [8] (see also [7],[15]-[16]) 



cos 2 (or) , 



V(4>) = - cosh 2 



Yf( c2 -# 2 ). 



0= L/3(c2_if2) ln 



3# 2 + c 2 -4c 2 tanh 2 



1 + tan 



tan 



V 



H 2 ) 



(38) 
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where c is an arbitrary constant. The 5d manifold with this choice of A(cr) has naked 
singularities at a — ±7r/(2c) and the scalar field diverges at this singular point. However, 
in the particular case c = H, the scalar field vanishes and the potential becomes constant, 
yielding the action (1) with the corresponding solution (7). 

It should be mentioned that when the author of [8] analyses the dynamics of metric 
fluctuations, he does not properly define the mass in a 4d de Sitter spacetime and, conse- 
quently, the corresponding Schrodinger equation possesses an analog quantum mechanical 
potential which vanishes asymptotically when c — H, however, he correctly reproduces the 
pair of bound states present in the model for this limit. 

On the other hand, Wang [15] has obtained a family of solutions for the same field 
system (37), which was borrowed from a 4d solution presented in [21], and in terms of the 
z coordinate reads: 

A = -n In [cosh (cz)} , <p = <po sin' 1 [tanh (cz)] , V{<j>) = V cos 2(1_n) (-£-), (39) 



where 0o = y/3n{l - n), V = nc 2 [3(l + 3n)]/2, c and n are arbitrary constants and the 
latter is subject to < n < 1. 

Because of this inequality, in principle, this family of solutions does not contain our 
particular solution as a special case, but as a limit n — >■ 1, when the scalar field vanishes 
and the potential becomes constant. 

The 5d curvature scalar for this solution reads 

R = 4c 2 n(3n + 2) cosh 2(n ~ 1} cz (40) 

which is regular for n < 1 and diverges for n > 1; it is worth noticing that it becomes 
constant in the limit n — > 1, a fact which agrees with the analysis performed in the previous 
section regarding the singular or regular character of the 5d manifold. 

Finally, it is worth noticing that corrections to Newton's law were not computed either 
in [8] nor in [15]. 
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